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Abstract 

In this paper we present the Koopman-von Neumann (KvN) formulation of classical non- 
Abelian gauge field theories. In particular we shall explore the functional (or classical 
path integral) counterpart of the KvN method. In the quantum path integral quantization 
of Yang-Mills theories concepts like gauge-fixing and Faddeev-Popov determinant appear 
in a quite natural way. We will prove that these same objects are needed also in this 
classical path integral formulation for Yang-Mills theories. We shall also explore the 
classical path integral counterpart of the BFV formalism and build all the associated 
universal and gauge charges. These last are quite different from the analog quantum ones 
and we shall show the relation between the two. This paper lays the foundation of this 
formalism which, due to the many auxiliary fields present, is rather heavy. Applications 
to specific topics outlined in the paper will appear in later publications. 



1 Introduction 



In 1931, soon after the appearance of quantum mechanics (QM) in its operatorial formulation, 
Koopman and von Neumann (KvN) proposed 1 an operatorial version also for classical 
mechanics (CM). It is well known that, given any operatorial formalism, it is always possible 
to build an associated path integral, like Feynman did in 1942 for QM. The same was done 
in 1989 T for the KvN operatorial formalism of CM. We will indicate this path integral with 
the acronym CPI for Classical Path Integral in order to distinguish it from the Feynman's 
one associated to QM which we will indicate with QPI for Quantum Path Integral. For a 
brief review of the CPI of a point particle we refer the interested reader to Ref. [2j and later 
papers 0. In this section we want to briefly review the extension of the CPI to the case of 
a system with an infinite number of degrees of freedom, i.e. the scalar field theory with a ip'^ 
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interaction already introduced in This case is free from the technical comphcations 

due to the constraints present in a field theory with gauge invariance. The Lagrangian density 
describing a scalar field (p with a 99^ self-interaction is given by: 

1/^ N/^,,^ 1 00 1 A 



79^ 



1.1 



Indicating with ir the conjugate momentum of f we can build from the Lagrangian Hl.l|l the 
following Hamiltonian: 



H 



dx 



2 2^ ^ 2 24 



where dx means the integration over all the space variables and dr'Ji = t; — stands for the 

ax*" 

space derivative of the field ip. The equations of motion of the system, which are 
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can be written as a couple of first order equations in terms of the Hamiltonian H as follows: 

:i.2 



Lp{x) 



dH 



7r{x) 



dH 



dTr{x) dip{x) 

where d indicates a functional derivative and x collects together all the space-time variables 
x = (t,x). Using the compact notation = (97, vr), a = 1,2 and the 2x2 symplectic matrix 
00°"^ we can write (11.21) as: 



ab 



dH 



b • 



:i.3) 



We shall now implement the standard steps of the CPI procedure Let us ask ourselves, if 
we start from an initial field at time tj, which is the probability amplitude of finding the 
system in at time t? This probability amplitude is one if the field lies on the solutions 
of the classical equations of motion and zero otherwise. To write this probability amplitude 
as a path integral we must integrate over all the paths in phase space with a weight which is 
given by a functional Dirac delta: 



ir(r(x),tier(x),i. 



(1.4) 



where ^^[(x; ^"(x), tj) is the solution of the equations of motion (|1.3|) with initial condition 
(,"'{ti,x.) = ^"(x). The functional Dirac delta on the solutions of the equations of motion can 
be converted into a functional Dirac delta on the associated equations, provided we introduce 
also a suitable functional determinant 



e - co^^dbH 



det 



5(x - y)6^,dt - 



ab 



d^H 



5e(x)ae(y) 



:i.5) 



If we use the integral representation of the Dirac delta we can rewrite the Dirac delta on the 
RHS of (II3J) as: 
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- uj^'^dbH 



VXa exp 



i I dx\a{x) e{x)-u;^''dbH{0 
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The functional determinant in ()1.5() can be exponentiated via suitable Grassmann fields c° (x) 
and Ca(x) using the formula 



det 



5(x - y)5,"5t - u 



ab 



i dxcaix) 5',dt-uj'''dbdcHiO c%x 



VcDc exp 



The final result of these formal manipulations is that the probability amplitude H1.4|l can be 
rewritten as the following path integral over an extended phase space given not only by the 
fields 4 but also by the auxiliary fields (A,c, c): 



i dx C 



K{C,t\Cf,ti) = J V^VWcVcexp 
where C = Xa^"" + ^Cqc" — TC and 7i is the following Hamiltonian density: 



:i.6) 



:i.7) 



From a formal point of view the kernel of propagation H1.6|) is identical to the one given in [2] 
for a point particle. The main differences are due to the fact that a field is a system with an 
infinite number of degrees of freedom. This implies that the derivatives Sf, must be intended 
as functional derivatives and that in the weight of the path integral the Lagrangian density 
C is integrated over all the space-time variables x. Another difference is in the functional 
measure: for example in the expression 2?^" = Hfc x '^^°(^' x)' besides a slicing in t, there is 
also a slicing in the space variables x. The CPI for the point particle has been thoroughly 
studied |21- 0-101 from a geometrical point of view and contact has been made both with the 
KvN formalism and the differential geometry approach to CM. In particular, Ti turns out to 
be the Lie derivative 7_ of the Hamiltonian flow associated to H. 

The reader may wonder why we want to do the CPI for a field theory. The reasons are 
basically two: the first one is that some cZassica/ solutions (solitons, instantons, monopoles, 
etc.) seem to play an important role in several non-perturbative aspects of quantum field 
theory, like confinement, phase transitions in the early universe and similar (9j. So a formalism 
like the CPI, which is classical but formally quite similar to the QPI, may help in unveiling 
further connections between classical solutions and quantum effects. The second reason for 
the CPI is that classical field theory seems to play a central role in heavy-ion collisions |inj . 
The reason is that in these collisions quantum effects are damped, due to the very high 
occupation numbers. People ^J-ffS have realized that even the perturbation theory to be 
used in heavy-ion collisions should be the classical and not the quantum one. Clearly the 
best tool to use in order to derive the Feynman diagrams for classical field theory is the CPI. 
This CPI perturbative approach was first pioneered in Ref. for a (p^ scalar field theory. 
Of course the classical fields involved in heavy-ion collisions are gauge and fermion fields. So 
a CPI version of these fields should be developed. For fermions things were done in Ref. 0, 
while for gauge fields we shall work things out in this paper, where we will concentrate on 
the more formal aspects of the CPI for YM. In particular, we will study the role that the 
gauge-fixing and the Faddeev-Popov determinant play even at the classical level both in the 
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Faddeev-Popov (FP) and in the Batalin-Pradkin-Vilkovisky (BFV) jl3j approach. Like in 
the standard CPI for the point particle we have, due to the many auxihary variables, a lot 
of symmetries (like super symmetry, ghost charge invariance etc.) which are universal in the 
sense that they are present for every Hamiltonian H entering (|1.7j) . In the CPI for Yang-Mills 
(YM) theories these symmetries will be even more numerous, due to the gauge invariance of 
the theory. These formal tools will be developed in this first paper. In a second one we 
will implement the perturbation theory both at the scalar field theory jl2j and at the YM 
level. In particular, thanks to the supersymmetry |2j that we will work out in this first paper, 
we shall show |14| that the complicated ^2] perturbation theory of the CPI can be simplified 
using the analog of the superfields studied in ^^1- This perturbation theory will turn out 
to be very similar to the quantum one because, once we replace fields with superfields, the 
QPI goes into the CPI [Kj . We hope this simplified perturbation theory can help the nuclear 
physicists \1'2\ who are exploring the quark-gluon phase in the heavy-ion collisions. 

The supersymmetry studied in this first paper will be the central tool which we shall use 
also in the third paper ^1 of this series which we hope to dedicate to the topological field 
theory aspects of the CPI for YM. It was in fact proved in [2j that any CPI can be turned 
into a topological field theory by properly choosing the boundary conditions. There it 
was shown that the "topological" phase of the CPI helps in building objects like the Euler 
number of the phase space or the Maslov index of the Lagrangian Grassmannian associated 
with the phase space. In the CPI for YM the analog of the phase space will be the space 
of gauge orbits, so the associated topological field theory will help in building the analog of 
the Euler number or Maslov index of the space of gauge orbits. These indices will help in 
unveiling some of the geometrical features of the space of gauge orbits which are known to 
be at the heart of many important non-perturbative aspects of gauge theories. 

Before concluding this Introduction, we should also remind the reader of the lot of work 
which has been done since the early 1970's jl7j on classical field theory. In those papers the 
classical tools of symplectic differential geometry were put at work for classical field theory. 
Those tools are the same which enter our CPI. It was in fact shown in that all the 

several auxiliary variables of the CPI have a nice geometrical meaning. For example, some of 
the Grassmann variables can be interpreted as differential forms, the charge associated to 
one of our universal symmetries turns out to be the exterior derivative and the Hamiltonian 
Tl associated to C is basically, as we already said, the Lie derivative |7j of the Hamiltonian 
flow. These common tools can help the reader familiar with jjj to easily understand the CPI 
P] and vice versa. 

The paper is organized as follows. In Sec. 2 we review the YM formalism from the 
constraints point of view and build the associated CPI. This first version of the CPI will be 
called the "natural" one because it arises in a very natural manner from simple considerations. 
In Sec. 3 we review the BFV method J3| and build the associated CPI. In Sec. 4 we prove 
the equivalence between the "natural" CPI introduced in Sec. 2 and the BFV one of Sec. 3. 
In Sec. 5 we derive all the universal symmetries present in any CPI but now adapted to YM. 
We also build the charges related to the gauge invariances which have a different expression 
than in the quantum case because of the many auxiliary fields present. In Sec. 6 we build the 
universal supersymmetry present in this formalism and put forward some ideas for further 
developments which will appear in Ref. 1161. Several calculational details are confined to 
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few appendices. Because of tlie nature of this Journal we have tried to present things in a 
review form reviewing some topics which are somehow scattered in the literature. We have 
also presented the calculations in details because this first paper should contain many of the 
technical tools which will be used in further development of this project. 

2 Yang-Mills Theories 

As we have seen in the previous section, the implementation of the CPI for a classical scalar 
field theory is straightforward. Things become much more difficult when we study YM 
theories because of the presence of constraints. First of all, let us review some aspects of 
YM theories, just to fix the notation we will use later on. YM theories generalize quantum 
electrodynamics to systems endowed with a non-Abelian gauge invariance. The generators 
of the symmetry transformations form a Lie group Q, so they obey a well-defined algebra. 
If with Ta, a = 1, . . . , n we indicate the basis elements of this Lie algebra, we will write the 
composition law among them as: 

where C^^ are the so called structure constants. The antisymmetry and the Jacobi identity 
of the commutator imply that also the structure constants are antisymmetric and satisfy the 
Jacobi identity: 

ah — ~'^ba i '-'ab ^ec "r ^fec ^ea "r ^ca ^eb — ^■ 

Next we can define the second rank tensor g^b = C^^C^^^, which can be used to raise and 
lower the color indices. Another convention that we must specify is related with the choice 
of the space-time metric. We will use the following one: 



/-I 
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to raise and lower the position of the space indices. For example, on the four vector and on 
the space-time derivatives the action of the metric tensor is the following: 

= {t, x) =^ Xf, = -q^^x" = {-t, x), 

The gauge transformation on the field A has the following form: 

SAiiix) = AT{x) - A^,{x) = C,'^ebix)A>^{x) + (2.1) 

where g is the coupling constant, which, from now on, we shall put equal to one. If we define 
the antisymmetric field strength: 

Fr^d'^A:-d''A^,-gC,''^A^A'^, 
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we can deduce immediately from Eq. H2.1|) that the gauge transformations of this field are 
5Fa^ = C^'^e^Fc'^ . The standard gauge invariant Lagrangian density is: 

r — _1 7?M!^/?a 
From C we can derive the following YM equations: 

which can be written in terms of the covariant derivative D^^^ = 5^dy — C\^Al, as follows: 

^ a.^r = 0- (2-2) 
It is interesting to rewrite Eq. (|2.2|) in components: the component /j, = becomes: 

dkE'^ - Ci^E^Al = 0, (2.3) 

where we have defined the analog of the electric field E^ = F^^. Eq. ()2.3() is the analog of 
the Gauss law for electromagnetism. Since they do not contain time derivatives, Eqs. (|2.3|) 
are just constraints^. The space components of Eq. 1)2.2(1 are instead real equations of motion 
and are given by: 

-^^a + V X + CiJbA', + Ci,Bh X 1'^ = 0, 

where B is the analog of the magnetic field and its components are defined as 5^ = -e'^^F^^. 

Like for the scalar field theory, in order to implement the CPI we need to build from 
the gauge invariant Lagrangian density C = —-^F^^F^^ the associated Hamiltonian. In this 
case we must take into account the constraints of the theory. According to the calculations 
presented in Appendix A, the canonical action is given by |13| : 

S, = y" dx • vf^ - ivr^ • vf^ - ^4 • - Al{-V ■ vf^ + C^vf^ A')^, (2.4) 

where vr^ = —F^^ and the arrows above indicates the space components of the fields. In 
the previous formula, all the scalar products involve only the space components of the fields. 
From ()2.4|) we identify immediately the following YM Hamiltonian density: 

and the constraints: 

<ya = -y-^a + C\,^b • A' « 0. (2.5) 

Note that the A'^ in (|2.4|) act as Lagrange multipliers which will be indicated later on with 
A". The constraints a form a closed algebra, i.e.: 

{o-a(t,x),cJb(t,x')} = C"=^bf^c(t,x)5(x-x'). 



^As we will prove in Appendix A, they are the equations of the secondary constraints. 
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The above equation tells us that the Poisson brackets among the constraints are zero on 
the constraint surface, which implies that, according to the definitions of Dirac ^Hl) the 
constraints are first class. With similar techniques it is easy to prove that the YM Hamiltonian 
has zero Poisson brackets with the constraints: 



|aa(t,x'),y" dxHO^(t,x)| =0. 



This means that the Ga are automatically conserved in time and no other constraint is gen- 
erated. 

From the canonical action (|2.4|) we can derive the equations of motion. By considering the 
variation with respect to A'^ we get: 

Hi = ^ = -d,F',^ - X'&.A + CiAiF,'\ (2.6) 

Similarly, considering the variation with respect to tt^, we get the equation of motion for A^: 

i^ = < + 5fcA'^ + ClA'^^^ (2.7) 
Now, if we identify from the canonical action the following Hamiltonian: 

H = y dxj ivf^ • vf" + ]^Ba ■ + A°(- V • + C\,^h ■ A') I 

we can rewrite the equations of motion 1)2. 6|) and 1)2. 7|) as: 

7n dH dH , , 

dTTa dA" 

Collecting together all the fields in the variables (/)^ = (A'^,TTa), and introducing the symplectic 
matrix uj'^^ ~ (^^^l o) ' equations of motion (|2.8j) can be written in a compact form as: 

where dg stands for the functional derivative w.r.t. 0^. 

Let us now proceed to implement the CPI. In order to do that we must have a one-to- 
one correspondence between the equations of motion and their solutions, once the initial 
conditions are given. This is clear from formula H1.5() . Its RHS is equivalent to the LHS only 
if, given some initial conditions, a unique solution can be found for the equations of motion 
appearing in the Dirac delta on the RHS of (|1.5|) . If the solutions or the "zeros" of the 
equations of motion of the RHS of (|1.5)) were more than one, one should use the functional 
analog of the following formula: 



dl 




dx 
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where the sum on the LHS is over all the zeros Xi of the function f{x). The sum present 

df 

on the LHS in (|2.9|) does not allow us to bring the determinant 



dx 



to the RHS, as it 



is needed in order to implement the CPI. A similar problem has already been faced in the 
case of the point particle with reparametrization invariance • In the case of a theory with 
gauge invariance, like YM theories described by the action (|2.4j) . the number of solutions of 
the equations of motion with the same initial conditions is more than one. In fact, since the 
Lagrangian multipliers A" = are undetermined, when we solve the equations of motion 
we get an entire family of solutions depending on the parameters A", even if we fix the initial 
conditions. This is the reason why we cannot write the path integral, using the same formal 
expression of the scalar field H1.5() . as: 



P0^5(0^-cj^^5s//)det 



(2.10) 

Before passing from the Dirac delta of the solutions to the one of the equations of motions, 
we need to determine the Lagrange multipliers. This is possible by fixing the gauge jl8j-|19j. 
Suppose we choose the Lorentz gauge —d^A*^ = 0. To implement it we can use the standard 
procedure [2()j which consists in multiplying Z™^ by "one" and in rewriting this "one" as: 



1 



VAo 6[-d''A^]A^p{A), 



where App{A) is the FP determinant associated with the gauge-fixing —d^A^j^ = 0. In this 
way, Eq. 1)2. lUj) becomes: 



(2.11) 



•det 



V(l)^VAo 6[-d^A''^]App[A] 6{^^ - w-^^Ssi^) 



d(x - y)6^dt - 



The appearance of the FP determinant in a classical context is not a surprise: in fact, 
it is possible to prove that such a determinant can be reduced, via a suitable canonical 
transformation, to the determinant appearing in the Dirac procedure for classical systems 
with constraints, see for example IS or 21.. 

Following , we want to find the explicit form for the FP determinant associated with the 
choice of the Lorentz gauge. First of all, according to (|2.1j) . the generic gauge transformation 
of the field A is given by: 

6A^,=C,'%A^ + d>^ea 

which, for a gauge transformation in a neighborhood of the identity, implies that the gauge- 
transformed —d^A^j^ is given by: 



From this expression we get that the FP determinant is given by: 
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This determinant can be exponentiated via a couple of Grassmann variables tpa^Va'- 
We can exponentiate also the gauge- fixing via an auxiliary variable^ iTa- 



ViTa exp 



Like in the standard CPI procedure, also in this case the Dirac delta of the equations of 
motion, which appears on the RHS of ()2.11|1 . can be Fourier transformed via some bosonic 
auxiliary variables A^. The associated functional determinant can be exponentiated via a 
couple of Grassmann variables T'* and T^. The final result is: 



7YM(NAT) 
^CPI 



P/i exp 



dx e 



where the functional measure P/x is given by: 
while the Lagrangian C^^'^ is: 



(2.12) 



(2.13) 



UJ 



■ UJ 



'dsdcH)T^. (2.14) 



What we have done up to now is the first, more natural procedure to derive a path integral 
for a classical field theory with a non-Abelian gauge invariance. This is the reason why we 
have put an acronym "NAT" (for Natural) everywhere in the previous formulae. 

We want to conclude this section by underlying again that, in order to implement the CPI, 
it was crucial to fix the gauge and as a consequence to introduce the FP determinant. So a 
gauge-fixing procedure is needed even at the classical level and not just in the quantum case. 
In quantum mechanics the gauge- fixing was needed in order to get the quantum propagator, or 
equivalently to be able to invert the operator in the kinetic part of the Lagrangian. Somehow 
also at the classical level we need the gauge-fixing for a similar reason, i.e. to "propagate" 
the initial conditions in a well defined manner, by choosing only one of the many different 
trajectories. Even the FP determinant enters the CPI for reasons identical to the quantum 
case. This is very clear in our formalism while in the usual Dirac |18j procedure, the FP 
determinant makes its appearance via the "C" matrix needed to build the Dirac brackets 
from the Poisson ones. 

In the next section, we will implement another path integral for classical YM theories, 
which consists in starting from the so called BFV Lagrangian and in exponentiating the 
associated classical equations of motion. Even if they seem to be completely different, the 
two CPIs, the one analyzed previously and the BFV one which we will study in the next 
section, are instead totally equivalent, as we shall show later. 

■^The TTa should not be confused with the yfa appearing in 12.41 which carry a space index tt^. 
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3 BFV Method 



In the FP method analyzed in the previous section we seem to lose the original gauge in- 
variance of the theory because of the gauge-fixing procedure. Nevertheless, the functional 
integral maintain a symmetry, which carries the memory of the original gauge invariance of 
the theory, and which is called BRST symmetry, see for example ^HI- This symmetry is 
crucial because it guarantees that nothing observable depends on which gauge-fixing we have 
used. At the end of the Seventies, an alternative approach to the FP procedure was developed 
by Batalin, Fradkin and Vilkovisky (BFV). Their method gives automatically a gauge- fixed 
action with all the extra terms of the "FP type" plus the BRST symmetry mentioned above 

In what follows we briefly review the BFV method, following Refs. [^Hj- The nice 

feature of this method is that it is entirely classical even if it is mostly used at the quantum 
level. It is classical in the sense that it basically produces in an automatized way a gauge- 
fixing (GF) and the extra terms of the FP type which are things that are needed even at the 
classical level, as shown in the previous section. The BFV method starts by associating a 
couple of auxiliary variables V and Va to each first class constraint da and by defining the 
Poisson brackets among them: 



All the Poisson brackets between the Grassmannian odd variables {V,V) and the original 
phase space variables are instead identically equal to zero. We will assign to the Grassmannian 
odd variables a grading factor, called ghost number: in particular V will be assigned ghost 
number while Va will be assigned ghost number —1. Furthermore, the Lagrangian 
multipliers A'* = which enter the definition of the canonical action (|2.4I) must be considered 
as dynamical variables, so we associate to each of them a bosonic conjugate momentum vTa, 
which satisfies the following Poisson brackets: 



{vra,A^} 



5l (3.1) 



Of course, if we do not want to change the dynamical content of the theory, these momenta vr^ 
must be zero and this happens if does not appear in the Lagrangian. To summarize, if we 
indicate with the canonical variables describing the original gauge fields, the full extended 
BFV space is given by {z^, X°',TTa,V ,Va)- It is possible to prove that in this enlarged BFV 
space there exists a real Grassmannian odd charge which satisfies the following properties: 

ghost number ofO = +1 (3-2) 
= 0. 

The charge which satisfies the previous properties is the so called BRST charge: 

17 = / dx Tda - \v''V"C%Va 
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In Dirac's theory of constraints the observables are those functions of the original phase space 
which have zero Poisson brackets with the constraints, once we sit on the constraint surface. 
This means that the Poisson brackets of an observable O with every constraint aa must be a 
linear combination of the constraints themselves: {O, aa} = Wj^ai, ~ 0. We call "observables" 
the O because they remain invariant under the gauge transformations generated by the first 
class constraints cJa- Furthermore, two observables O and O' are equivalent if their difference 
is just an arbitrary linear combination of the constraints da, i.e. O' = O + k'^^a- Let us see 
which is the analog of these definitions in the BFV method. First of all, a function is an 
observable if it commutes with the BRST charge. The procedure to construct the observables 
is to extend the functions O to functions O of the BFV space by requiring that: 

• the functions O reduce to O when the ghosts V'^ and the antighosts Va are put equal 
to zero; 

• the functions O have ghost number zero; 

• O are real and Grassmannian even. 

• O are BRST invariant. 

Now, if O is an extension of the function O which commutes with the BRST charge and sat- 
isfies the previous properties, it is possible to consider an entire class of equivalent extensions 
Oe which differ from O by a pure BRST variation: Og = O — {9,Q}. Since the Hamiltonian 
itself is a particular observable, it also admits an entire family of BRST invariant extensions: 

He = H - {9,n}. (3.3) 

The function 9 is somehow related to the gauge freedom of the theory: fixing the gauge means 
choosing a particular function 9 and consequently a particular Hamiltonian among the Hg 
of Eq. ()3.3() . Clearly, all the physical quantities must be independent of the choice of this 
function 9 and this is guaranteed by the so called Fradkin and Vilkovisky theorem. 

Let us now apply the BFV method to the case of YM theories. As we already said before, 
the Lagrange multipliers A" can be identified with the time components of the fields A^, i.e., 
A^, so the associated conjugate momenta must be identified with the time components of the 
momenta n^. The equal-time Poisson brackets are: 

{7T,{t, x), \-{t, y)} = -5,"<5(x - y). (3.4) 

From the canonical action ()2.4() we get the following constraints vTa ~ 0, a = 1, . . . n. We 
know from Appendix A that they are the primary constraints of the theory. From their 
consistency, i.e. the requirement that their time evolution leaves them invariant, we get just 
the n secondary constraints of Eq. ()2.5() . so the total number of constraints is 2n. We will 
indicate them as 

o-fl = (vrfe,cjfe), a = l,...,2n; b=l,...n. (3.5) 

Following the BFV procedure, we shall associate to these constraints the ghosts Va and P". 
All these auxiliary variables are Grassmannian odd and, as we have seen, they are conjugated 
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in the sense that: {Pa^V } = -5a- These ghosts can be divided in two blocks according to 
the type of constraints H3.5|) they refer to: 

T = {-iP\&), ra = {iCb,Pb) 

and their Poisson brackets are given by: 

{P,,C''} = -6'i = {P\Ca}. 

Since the only non-zero structure constants are the ones among the secondary constraints, 
the BRST charge for the Yang-Mills theories becomes: 



dx 
dx 



UaC - iP^T^a + \PaCl,C''C' 



(3.6) 



The gauge freedom of the theory is entirely contained in the function 6 which appears in 
(|3.H)) . Usually this arbitrariness is moved to a new function x according to the definition: 

This function x plays the role of the GF jl3j . For example, the analog of the Coulomb gauge 
can be obtained by choosing x° = V • A". With all these ingredients in our hands, we can 
now proceed to build the BFV action in the form given in Ref. jl3j : 



dx 



where if 6 is built out of the x function associated with the "Coulomb-like" gauge mentioned 
above then the commutator {9, $7} is given by: 



dx 



-TTaV ■ A'' - X^aa - iPaP" 

With this choice the specific form of S^^^ is: 

5^^^ = j dx[ia ■ Vf^ + A'^vr, + P'^Ca + C'^Pa " ^vf, • Tt" - ^F'J F.^^ - TTaV ■ A'' + 

-X'^C^ac^b ■ A' - iPaP" + X'^PbCi.C + iCaV ■ (VC7« + ClJ'C'')}. 



(3.7) 



(3.8) 



The previous expression is BRST invariant and from it we can read out the following BFV 
Hamiltonian: 



1 



1 



dx <; ^Tfa • + iF'jF,^^ + 7r„V • A^ - A'^V • vr^ + A^C^vrfe • A 



iPaP" - X^PbCi^C - iCaV ■ (VC + Cl^A^C 



(3.9) 
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From the BFV action p.Hj) we can derive, by variation with respect to all the fields of the 
theory, the following equations of motion: 



A" 
A" 

pa 

'Pa 

^a 



-(V X 



vf^ + VA" + C'IiXA'' 



Ci,{Bb X A^) + VTTa 
1 1 



iPa 

iV[(W,"-Cr^,)C7^] 

-\''PcC\, - iV ■ VCa + iCUvCd) ■ A'' 



(3.10) 



/^a \b/^c 



If we introduce a compact notation indicating with all the fields above 

e^{A\\^^a,7Ta.P\C\Ca,Pa), 

then Eqs. (|3.1()j) can be written in the following form^: 



(3.11) 



where the index A in runs over all the fields of the theory, and we have used the symplectic 
matrix iv^^ , which has the following expression: 





1" 


A" 






pa 


(ja 


Ca 


Pa 


A" 


{ 





1 














\ 













1 
















-1 























TTa 





-1 




















pa 




















-1 





Qa 























-1 


Ca 














-1 











Pa 


I 














-1 





) 



The matrix above has an antisymmetric block associated with the commuting fields, and a 
symmetric one associated with the BFV Grassmannian odd fields. So the symmetry proper- 
ties of the matrix w are: u)^^ = —(—)''*' '^'a;^^, where \A\ indicates the Grassmannian parity 
of the field associated to 

Indicating with ujab the inverse matrix of uo^^ we can write in a compact way also the 
kinetic part of the BFV action. In fact, up to surface terms, we get: 



1 



1 

2 

A„ 



vf^ + A'^vr, + P'^Ca + C^Pa. 



^In the equation of motion 13.111 , which involves also Grassmannian odd fields, it is important to distinguish 
between right d and left d derivatives. If we want to pass from right to left derivatives and vice versa, we 



have to use the following formula: d bF — {- 
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Using the previous formula, the BFV action becomes: 



As we said previously, and as it is clear from Ref. [13], this BFV procedure is a classical 
construction (and not necessarily a quantum one), which has automatically generated the GF 
and the analog of the FP terms. Because it is classical, we can think of using its equations 
of motion in the construction of the CPI. In this case, as the BFV action contains already 
the GF, we can explicitly solve the equations of motion without encountering the problem of 
the previous section. This means that we get one solution out of a definite initial condition, 
so the CPI procedure can go through without having to add the GF and FP in the measure 
as we had to do in the Zcj^'^*^' of the previous section. 

So let us implement the CPI associated with the BFV action S^^^ . As usual, the initial 
expression is given by a functional Dirac delta over the classical solutions of the equations of 
motion: 



YM(BFV) 
CPI 



(3.12) 



In the next step we have to pass from 1)3.12(1 to the Dirac deltas containing the equations of 
motion. In doing that, we must take into account that we have both Grassmannian even [x) 
and Grassmann odd variables {9) among the fields Let us consider a generic change of 
variables 



A 
C 



dx 

where A = — — 



B 
D 

de 



dx do ^ ^ 

„ , , B = — — , C = — — - , D = — — . The blocks A and D are Grassmannian 
ox' oB' ox' od' 

even, while the blocks B and C are Grassmannian odd. The result of this change of variables 

in the integration measure is the following 



J] dx, n = / n • n k ^^^t 



A B 
C D 



(3.13) 



The expression sdet means "superdeterminant" |24j and its definition is given by 

sdet ^ = det{A - B D'^ C)det-\D) . 

Let us notice that, when only Grassmannian even variables are present, the superdeterminant 
reduces to an ordinary determinant while, in presence of only Grassmannian odd variables, 
the superdeterminant becomes, as usual, the inverse of a determinant. In the case of our 
CPI we can say that, in going from the Dirac delta of the solutions to the Dirac delta of the 
equations of motion, a functional superdeterminant appears: 



■ ijj 



a s^^'''' (6 sdet 



AC 



d 



. d 



(3.14) 
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The next step is to exponentiate the equations of motion and the superdeterminant. The 
equations of motion can be exponentiated by using the following formula: 



VK exp 



(3.15) 



where A is a field with the same Grassmann parity of the field it refers to: [A^] = [^^] = [A\. 
Also the superdeterminant can be exponentiated by using a generalized Gaussian integral, 
see Appendix B for details: 



sdet 



VTaT>T^ exp 



(3.16) 



The auxiliary fields T and F have a Grassmann parity which is opposite w.r.t. the one of the 
fields they refer to: [T^] = [f^] = [^^] + 1 = [^] + 1. 

Collecting all the terms in (HTT^ . (Hnil) . and (ITTHl) . we get: 



^YM{BFV) 



VA^V^^'VT^VT^ exp 



i / dxC^ 



where the Lagrangian density is given by: 



A. 



(3.17) 



(3.18) 



Since the previous Lagrangian density contains only first order time derivatives, we can 
immediately obtain the following Hamiltonian density: 



From C^^^ of 1)3. 18() we can derive the following equation for the fields F: 



dt6i-co^''dcHdB]r'' = 0. 



(3.19) 



(3.20) 



Like for the point particle P], also in the case of YM theories (|3.2Uj) is the equation of evolution 
of the Jacobi fields, i.e. the equation of evolution of the first variations of the theory, as we 
will prove in Appendix C. 

The BFV Lagrangian density and Hamiltonian density TC^^^ assume a particularly 
simple form if we combine all the fields of our theory in a set of real superfields: 



(3.21) 



where 9 and 9 are two anticommuting parameters and (— )''^' indicates the Grassmannian 
parity of the field The fact that such Grassmannian parity appears means that, differently 
from the point particle case, we have to distinguish between the even superfields defined as: 



-A 

■^even 



t + 9T^ + 9uj^^rB + muj^^As 
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and the odd superfields defined as: 



Let us notice tliat each S'* has the same Grassmann parity of the field it is associated 
with. Using the definition (|3.21j) we can derive (up to surface terms for the Lagrangian) the 
fohowing relations: 

J d^n''^^ = i J ded6H''^^{E), j d-^C"^^ = i j dodev^^i^). (3.22) 

The derivation of these formulae is not presented here in details because it is the same, 
modulo some grading factors, as the one contained in Ref. ,15,^ for the point particle. So, 
in order to pass from the BFV Hamiltonian i^^pv (|3.9|) and Lagrangian L^^^ of Eq. 

(TUni to W^^^ and respectively, it is sufficient to replace in H^''^ and L^^^ the fields 
with the superfields H"* and then perform the Grassmann integration i J dOdO. 

Since we will make an intensive use of these superfields later on, it can be useful to write 
their explicit form using the following symbols: 



At 


= AI + er^fe 




+ ioek^u 




= vr^' + er^* 




k 


A" 


= A" + er^" 


+ er^^ + mA^^ 




= vr, + er^^ . 


- ^Ta" - 


- ieeAxa 






+ 9Tp 

J^^ a 


+ mA-p 




= Ca- er^" 


+ 6T pa 


+ i69Apa 


pa 




+ ^c. 


+ iOOAj, 

O a 


fa 


= Pa- 9r^^ 




+ i69Ac^ 



(3.23) 



4 Equivalence of the CPIs for Yang-Mills theories 

In the two previous sections we have implemented two different CPIs for YM theories, the 
"natural" and the BFV one. The presence of more than one path integral both at the 
classical and at the quantum level is due to the fact that in field theories there exist different 
procedures to fix the gauge, like the BFV and the FP ones. In this section we want to 
convince the reader that the two different functional integrals that we have obtained for the 
CPI are equivalent. 

First of all, let us prove that the equations of motion that can be derived from the standard 
FP effective action^ and from the BFV action are basically the same. If we choose the Lorentz 
gauge —5'^^^ = 0, the FP effective action becomes f^ : 

= - ^ad^AI - id^^ad,.rf + id^^aC\jfA\, (4.1) 

*We call "FP effective action" the full action, with GF terms and FP ghosts that one gets in the quantum 
path integral if one uses the FP procedure |18 | .|22 | . 
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where tTq are the auxihary variables used to exponentiate the GF, ipa and ri"" are the Grass- 
mannian odd variables used to exponentiate the FP determinant and are the structure 
constants of the gauge group under which —\F'^^Fa^ is invariant. In the previous section we 
have written the BFV equations of motion in a symplectic form and, if we want to compare 
the two procedures, i.e. the FP and the BFV one, we should write also the equations of 
motion that can be derived from 1)4. 1|1 in a symplectic form. To do this, we must derive from 
1)4. 1|) the associated Hamiltonian. We can start from the definition of the energy-momentum 
tensor containing both Grassmannian even and odd variables: 

nrFP f)T^^ f)T^^ 

d{d^Al) d{d^ipa) d{d^r]^) 



In particular, the 00 component gives: 



1 

4' 



+7r,V • A'' + iViJa ■ - i^i^a ■ A'Cl,r^'. (4.2) 

From the definition of the conjugate momenta = — ■■ — we get the following equations of 
motion for A: 

Aa = -VAl + C^'^AkoA, + 7?,. (4.3) 
Similarly, the momentum conjugated to A^^ is given by: 

-a = (4.4) 



Furthermore, by construction: 



{f^uFH'' = l^^F^ - ■ r. (4.5) 



For what concerns the momenta conjugated to the Grassmann variables, we have: 

from which we easily derive the equations of motion for rf and '\\)a- 

if = -m^^ + C%^A\rf, V-a = iTT^-. (4.6) 

From the previous equations we note that the FP ghosts evolve just as the BFV ghosts: 

= -ipa + C1,X''C', ta = iPa, 

SO we can make the following identifications: 

C^T], C^V- (4.7) 
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Consequently, the variables P and P, which appear in the BFV formalism, can be identified 
with the momenta conjugated to the FP ghosts: 



P ^ vr^, 

Replacing (|0|) - (|^ into we have: 



P ^ TTrj. 



^ - ij u 2 



Up to surface terms, we can identify the previous expression with the following FP Hamilto- 
nian density: 



which coincides with the BFV Hamiltonian of Eq. if we use the identifications (|4.7|) - 

H4.8p . This immediately implies that from the action 

S = y dx la • Vr" + X^TTa + i^aTT^^ + ^"t^t,- - 

we can derive, by variation with respect to the fields, the same equations of motion p.lUf) 
that we got via the BFV method. So we can say that the BFV method is somehow a 
canonical version of the FP method. In fact, either starting from the BFV or the effective FP 
Lagrangian, we get the same equations of motion in the symplectic form. As a consequence 
from their "exponentiation" we will get the same weight for the CPI, i.e., = C^^ . Note 
that C^^ is not the weight that appears in Zq^i^^'^^ of Eq. (|2.12j) because this last CPI was 
obtained not from the FP effective action but from the standard YM action. 
Finally, we want to prove that the CPI implemented using the BFV procedure 



7YM(BFV) 
^CPI 



y "D/i' exp 



dx 



(4.9) 



and the one obtained using the "natural procedure" Z^j^'^*^' are equivalent. First of all, we 
shall prove in Appendix D that 



(4.10) 



In this proof we made use of the equations of motion because in the CPI the system is always 
"sitting" on the classical trajectories. Eq. (|4.1()|) can help us in giving a further proof of 
£BFv _ £FP^ fact, using the superfields H of Eq. 1)3. 21() . we can pass from the usual BFV 
Lagrangian density L^^^ to the one which appears in the weight of the BFV classical path 
integral C^^^ according to the following equation: 



(4.11) 
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The same holds for the FP effective action: 

C^'' = i [ d0d0L^^(H). (4.12) 



So from L''^^ = L^^^ and Eqs. KTi\ and (BT2I1 we get again that £^^^ = ^C^^. This £^ 
appears in the weight of Zq^^ ^^^: 



■^cpi'^''' = J Vfi' exp i J dx C 



(4.13) 



where P/i' means the functional integral over all the fields present in the Lagrangian density 
C^^ . Since = C^^ and since the fields in the measures can be identified, as done in 1)4. 7() . 
(|4.8|) . we can say that Zq^^^^^^ = Z^^i^^\ Because of this, if we now want to show that the 
BFV classical path integral Zq^^^^^^ is equivalent to the ^^naturaF one Zq,^''"'*^' it will be 
sufficient to prove that the functional integral Zq^i^^^ in (|4.9j) is equivalent to the one given 

byEq. frm -. 

'dxC''^'^ . (4.14) 



^YM(NAT) I 

Zcpi ' = Vfi exp 

The comparison is easier if we divide the weight of the FP classical path integral in three 
parts: the GF, the ghost and the bosonic part: 

T^FP T^FP I T^FP I 7-FP 

The three parts can be written in terms of the superfields as: 

C^J = ij dOde {-nad'^A'l^} 

£-11 = i j dode [-id^^ta&^c,,] (4.15) 
cz, = I j dode i^-\F-^,F^^Y 

As usual, we have indicated with a tilde sign the superfields of Eq. (|3.23)) . Similarly, we can 
divide the natural Lagrangian in the following three parts: 

The main obstacle in proving that the two functional integrals (|4.13|) and 1)4. 14|) are equivalent 
is that the number of variables in C^^ is larger than the number of variables in in 
particular for what concerns the GF and the ghost parts. For example, the generating 
functional Z^^^^^^'' lacks all the functional integrals over the CPI ghosts F, F and A associated 
with the BFV ghosts. However it is possible to prove that all the extra fields present in Z™^^^^ 
can be integrated away to get just Zq^i^^'^\ Since the proof is rather long, we have confined 
it in Appendix E. 

So we can conclude this section by saying that the two different CPIs for YM theories that 
we have implemented in the previous sections, the BFV and the ''^naturaF one, are equivalent. 



19 



5 Universal symmetries of the CPI for Yang-Mills theories 



In the previous section, we have proved that the CPIs derived using the "naturar and the 
BFV procedures are equivalent, so in the study of the symmetries of the YM CPI we can 
concentrate ourselves on ()3.17() . which was obtained by implementing the BFV method. First 
of all, let us derive from (|.S.17|) the graded commutators of the theory. To do that we use 
the standard method developed by Feynman P3| for the standard commutators of quantum 
mechanics. In our case from ()3.17p we get: 

[e^(t,x),A^(t,y)] =<5(x-y), [r^(t, x), r^(t, y)] = 5f <5(x - y). 

All the other commutators are zero. In particular, the commutator among the fields ^ is zero: 
C^] = 0- We know that the CPI involves a number of fields which is larger than the one 
needed to describe the system, so, as it emerges from this redundancy in the description 
of the system implies the presence of a certain number of symmetries. We can immediately 
generalize those symmetries from the CPI of the point particle to the one for YM theories. 
For example, the BRS charge, derived in [2] for the point particle, becomes: 

Q{t) =i j dxr'*(t,x)A,4(t,x). (5.1) 

As usual, Q is Grassmannian odd and nilpotent and this is the reason why we called it a BRS 
charge. In spite of its name, it has nothing to do with the BRS charges |2(ij of gauge theories. 
In fact the BRS charge of the CPI made its appearance even in the point particle case j2] 
where no gauge field was present. There [2] we showed that the CPI-BRS charge could be 
identified with the exterior derivative [7j on the phase space. In the case of YM theories the 
phase space is built out of the whole variables (which include the gauge ghosts besides 
the gauge fields), so the BRS charge 1)5. 1|) becomes the exterior derivative on this enlarged 
space^. In the point particle case there was another Grassmannian odd and nilpotent charge 
which was called F21 the anti-BRS charge and that in the CPI for YM theories becomes: 

Q{t) = -i j dxA^(t,x)u;-"^rs(t,x). (5.2) 

The charges Q and Q anticommute among themselves and commute with the Hamiltonian 
/ dx'H^''^, i.e. they are conserved charges. 

In the CPI for a point particle there were three further conserved charges K, Qf, which 
were bilinear in the ghosts and that together with Q and Q made an ISp(2) algebra. What 
can we say in the case of YM theories? The expression P] of those charges must be slightly 
modified in the gauge case in order to preserve the ISp(2) algebra. They are: 

K = - J dK Ir^wsAP^, K = J djc ^r^o^^^r^, Qi = - J dxr^r^. (5.3) 

These three charges, modulo central charges, close on the Sp(2) algebra: 

[Qi,K] = 2K, [Qi,K] = -2K, [K,K] = Qf (5.4) 

^We will later show that in the case of the CPI for YM we have, besides this universal BRS charge, another 
one related to gauge theories (see Sec. 5.1). 
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and they are all conserved. Furthermore, then' graded commutators with the BRS and anti- 
BRS charges reproduce the inhomogeneous part of the ISp(2) algebra: 



[QuQ] = Q, [Qu Q] = -Q^JK, Q] = 0, 

[K,Q]=Q, \K,Q]=Q. [^,Q]=0. (5.5) 

Like for the point particle, also in the YM case the symmetry transformations generated by 
the previous charges can be represented in a very compact way by introducing the superfield 
formalism. As we have seen in Eq. (|3.21l) . the superfield associated with is given by: 

= + (-)'^ier'* + (-)'-"i^a;^^rs + i(-)[^'^0a;^^As. (5.6) 

Besides connecting the weights of the classical and quantum path integrals according to Eq. 
(|3.22|) . the definition of the superfields (|5.6)) allows us to represent the generators of the 
ISp(2) algebra in terms of differential operators ,2^ which act on the superspace (t, 9, 9) in 
the following way: 

The previous operators generate the correct transformations according to the standard rule: 

where VL is one of the operators of Eq. ()5.7|1 . These operators satisfy the same algebra of the 
associated charges. Before going on, let us notice that the operatorial representation of the 
ISp(2) charges on the superspace in the YM case is the same as the one of the original point 
particle formalism [^j, provided we change the definition of the superfields by introducing the 
suitable grading factors of Eq. (|5.6j) . 



5.1 The BRS-BFV charge 

In the previous sections we have used the superfield formalism to derive in a simple and 
elegant way TC^^^ and C^^^ from //^fv g^j^^ L^^^, see Eq. ()3.22|) . and to write in a compact 
way an operatorial representation for the generators of the ISp(2) algebra, see Eq. ()5.7|) . The 
superfields have also a third useful property: they allow us to write immediately a second 
conserved BRS charge. This charge is very different from the universal charges (|5.1j) - (|5.3|) 
that we have considered up to now. It is in fact related to the original gauge invariance of 
the theory and it exists only for the CPI of gauge theories. 

As we have already mentioned, even after the gauge-fixing procedure, in the quantum field 
theory a memory of the original gauge invariance still remains under the form of a symmetry 
generated by a gauge BRS charge 26 . This invariance is present at the QM level to guarantee 
that the theory is independent of the particular choice of the gauge-fixing we may have made. 
Something similar should happen also at the CPI level because also there we have a GF but 
at the same time we want everything physical be independent of which form of the GF we 
choose. To derive the form of the gauge BRS charge in the CPI formalism, let us remember 
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Eq. I|3.22() : there we have constructed the BFV-CPI Hamiltonian^ TL out of H by replacing 
the fields with the superfields and performing the integration over the Grassmann partners of 
time 6 and 9. Applying the same procedure to the charge of Eq. H3.6() . we get the following 
charge in the enlarged space of the CPI: 



With a calculation identical to the one performed for the Hamiltonian, we get: 

n = A^uj^'^'dB^ - iT^uj^^'dc^'dBT''. (5.8) 

Let us notice that, due to the Grassmannian odd character of Q, in (|5.8j> there appears a 
minus relative sign in the second piece. 

Now, let us consider two different generic functions, Oi and O2, of the standard BFV 
phase space with Poisson brackets'' {Oi, O^}^^^ = O3. As we will prove in Appendix F, the 
associated functions in the enlarged space of the CPI 15 : 

= A^c^-^^^sOi + iP^c^^^^flOi 9"cr^(-l)[°'l (5.9) 

satisfy the commutator 

[d„d,]=id,. (5.10) 

These functions Oi are nothing else than the Lie derivative of the Hamiltonian flow 

generated by Oi. The procedure to pass from Oi to Oi is analog to the momentum map |Ej 
or its inverse. As a particular case, let us come back to the BRS-BFV charge of Eq. p.6() . 
It is conserved in the standard BFV space because {0, H} = 0. From (|5.in|) this immediately 
implies that the commutator between the associated Lie derivatives is zero: 

[n,n] = 0. 

Since TC generates the time evolution in the extended CPI space, we can conclude that Q is 
a conserved charge. With similar calculations, the nilpotency of 0, implies the nilpotency of 

n-. 



= 



0. 



0, will be indicated as the BRS-BFV charge, to distinguish it from the other BRS charges, 
like Q, Q, present in the CPI of every model. We would like to stress that, without using the 
superfield formalism, it would have been very difficult to find complicated charges like (|5.8|) 
and to understand that they are symmetries for the Hamiltonian 7i. 

It is a long but straightforward calculation to prove that the graded commutator of 0, with 
the ISp(2) charges Q, Q, K, K, Q{ is zero. In particular, Q commutes with 17. This is crucial 
because [2j Q can be interpreted as an exterior derivative and we know that the exterior 

''From now on, we shall not indicate the dependence of the fields on x and the integral over x. We will 
also omit the index BFV on H and Ti.. 

d d 

^The Poisson brackets are defined as; {OijOg} 1, ~ Oitt— -o;'*^— — O2. 
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derivative is coordinate free, so it must be also gauge invariant. This is guaranteed just by 



the commutator 



0. The main steps of the proof are given in Appendix G. 

Before going on, we want to study the different physical meaning of the two symmetry 
charges Q and O. If we write in a table all the fields of the CPI obtained implementing the 
BFV method we get: 



Al TT^ A" 7r„ C« Ca Pa Pa 

r^fc r''" r^" r'^" r"^" pCapPa pPa 

^_Ai r^fe TTTa T_c- Tc'^fp^ fp^ 
~^Ai A^fe Aa„ AcKp^ Ap^ 



The action of the extended BRS-BFV charge Q on the fields of the first row via the CPI 
commutators is identical to the action of the charge via the BFV Poisson brackets. In fact, 
they both generate the gauge-BRS transformations which mix only the fields which appear 
in the first row of the table. The other BRS charge Q instead turns the fields of the first row 
into the associated Jacobi fields F, which are placed in the second row of the table, and the 
fields of the third row into those of the 4th. The anti-BRS Q instead turns the variables of 
the first row into those of the third row and the variables of the second row into those of the 
fourth. So if the charge allows us to move horizontally throughout the table, the charge 
Q moves us only in the vertical direction. We want to stress again that, even if we indicated 
both charges as BRS charges, the physical meaning of and Q is completely different: one 
comes from the gauge invariance, while the other comes from the pure CPI construction. In 
the original BFV theory |2Zj there were other conserved charges besides fi. By "original BFV 
theory" we mean the phase space and the Hamiltonian ff^pv ^[^\^ j^^g associated Poisson 
brackets. Those charges were an anti BRS-BFV charge: 

n = iPaTTa + aaCa + ^CaCbCJ^'P'' - PbC^C'Ca + TT, C^'^'^ A^'Ce , (5.11) 

a ghost BFV charge 

17g = iP^'Ca + iCPa, (5.12) 
and the analog of K and K which are given by: 

lC = -iP''C^ + ]^C,f\dC'Cf, lC = -iPdCd-\c'f''\dCfCe. (5.13) 

The five charges of Eqs. (|3.()j) . (|5.1H) - (|5.13|) satisfy another ISp(2) algebra. Out of these five 
charges we can build the associated CPI charges defined via Eq. ()5.9|) . let us call them il, fi, 
ilg, /C and /C. It is clear from Eq. (|5.1fl)) that these five charges close on an ISp(2) algebra 
in the enlarged space of the CPI and that they are conserved under the 7i of Eq. (|3.19|) . In 
Ref. ^ni; in connection with the topological field theory, we will also work out the various 
cohomologies associated to our Q and O charges. 
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6 Supersymmetry in classical mechanics 



In the second paper of Ref. H] it was proved that in the CPI for the point particle, besides 
the universal symmetries we presented so far, there was also a = 2 supersymmetry |28j . 
The charges associated to this supersymmetry depend explicitly on differently than the 
ones in Eq. (|5.5j) . We want to prove that it is possible to find the analog of these charges for 
YM theory provided we introduce suitable grading factors. First of all, let us notice that the 
following two quantities: 

are conserved because they commute with 7i. 

Let us now build the following linear combinations of N and N with the BRS and anti-BRS 
charges Q and Q of Eqs. (|^ - (f^ : 

Qn = Q-f3N, Q^ = Q + f3N, (6.1) 

where (3 is an arbitrary parameter. Being linear combination of conserved charges, Qh and 
Qu turn out to be also conserved. If (3 is infinitesimal and can be written in terms 
of Q and Q as: 

What is more important is that they satisfy the following anticommutator [Qh, Qh] = 2ipTi, 
so they can be considered as generators of non-relativistic supersymmetry transformations. 
Once we apply them on the Lagrangian, they generate the following transformations: 



and, as the Lagrangian changes by a total derivative, we can say that Qh and Q^ generate 
symmetries for the action S = J dt C. Finally, it is possible to prove that Qh and Qj^ are 
nilpotent: [Qh,Qh] = [Qh^Qh] = 0- 

As we did before with the BRS charges, we can give to these charges a differential operator 
representation in the superspace. Their action is defined via the equation 63-^ = —€0.3^, 
where is the YM superfield of Eq. (HT^ . The result is pH] : 

d-d — d d 

^ ^ d 

With the previous choice we have the following anticommutator [Qh,Qh] — 
— d 

operators Qh and Q,, close on — , which is the operatorial representation of the Hamiltonian 

at 

in the superspace. 

Every supersymmetric field theory is characterized by the number N of the independent 
supersymmetry charges whose square reproduces the Hamiltonian 

Qj = n, i = i,---,N. 
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We want to prove that in the CPI such number is given by = 2 j28j . First of all, let us 
consider again the four charges 

Ar = r^'a 4i/(_)M, N = TAUJ^''~d bH. 

Using their algebra: 

[Q,N] = [Q, N] = 0, [Q, N] = in, [Q, N] = -iU 

we can easily construct the "square roots" of the Hamiltonian 7i, which must have one of the 
following forms: 

Qa,/3 = aQ + I3N, J = aQ + pN, 

where a P, a and /3 are arbitrary parameters, which will be fixed in a while. Qa,i3 and Qa^p 
anticommute and satisfy: 

[Qa,p,Qa,p] = 2aPin, \Qa,-p^Qa,$]=-'iapm. (6.2) 

On the RHSs of (|6.2|) we get just the Hamiltonian 7i only if a/3 = — | and a(5 = |. Each of 
the previous equations can fix only one of the two parameters, so we have infinite possible 
square roots of the Hamiltonian. If we require that the charges commute with each other, we 
can extract the following four operators in the superspace: 

They are all conserved since they commute with the Hamiltonian. Furthermore, they com- 
mute among themselves and they can be considered as square roots of the Hamiltonian, in 

the sense that [Qj, Qi] = 2Q^ = Ti.. This implies that our supersymmetry is a = 2 super- 
symmetry, where Qi and Q2 can be identified with the supercharges and Di and D2 with 
the covariant derivatives jZW]. 

It is also possible to find the covariant derivatives associated with the original supersym- 
metry charges Qh ^ Qh- In analogy with what we have seen before, we can define: 

d-d — d d 

The previous operators commute with the supersymmetry operators Qh e Q^, so they can 
be considered as covariant derivatives. 

In supersymmetric field theories it is possible to associate with each covariant derivative 
a chiral superfield which is the one annihilated by the covariant derivative. Also in our case 
we can define a chiral superfield such that: Dn^ch = 0- From the previous equation we 
immediately get: 

EX = ^^ + {-)^^^u:^''WB-Oepi^. (6.3) 
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An antichiral superfield H^j^ is instead defined by the following equation: DfjE'^^ = 0, which 
is satisfied by 

The chiral and the antichiral superfield depend only on two types of fields. If we require 
that a superfield is both chiral and antichiral we get: = i.e., we reduce ourselves to 
consider only the fields of the standard BFV phase space. 

It is easy to show that if we apply a supersymmetry transformation on a chiral superfield, 
we get another chiral superfield. In fact let us indicate with Qi a supersymmetry charge, 
with Di the covariant derivative and with the associated chiral superfield. Then under a 
supersymmetry transformation the chiral superfield changes as follows: 

■=■011 = (1 ~ ^Qi) ^ch^ 

with e anticommuting parameter. If we apply Di to the transformed chiral superfield and we 
use the fact that [Qi,Di] = 0, we get: 

i.e. also is a chiral superfield. This means the space of chiral superfields carry a repre- 
sentation of the supersymmetry. Chiral superfields are important because they are made of 
fewer components than the standard superfields. This means that they carry a representation 
of supersymmetry of lower dimension than the one carried by the superfields. Moreover, it 
is possible to prove that their representation is an irreducible one |3Uj . 

The supersymmetry charges, as most of the elements present in the CPI, have a nice 
geometrical interpretation. It was proven in [2H1 that, like the BRS is an exterior derivative 
on phase space, the supersymmetry charges are equivariant exterior derivatives. We will 
not dwell here on this issue which has been thoroughly explored in Ref. |2BJ. The same 
geometrical meaning is for sure valid in the gauge case and this issue will be further explored 
in Ref. JH]) in connection to topological field theory. 

The supersymmetry that we have found in the CPI of any field theory is a universal 
one. It strongly indicates that from any non-supersymmetric field theory one can build, via 
the CPI, a supersymmetric one. So somehow supersymmetry seems to be a feature not of 
particular systems but of any system. Unfortunately this supersymmetry is non-relativistic. 
Why is so? We feel that this is due to the fact that, to build the CPI, we started from the 
Hamiltonian version of the equations of motion, which is intrinsecally non-covariant. One 
idea we would like to explore in is to start from that covariant version of the equations of 
motion, known as the de Donder-Weyl formalism, for a review see PJl. We feel that the CPI 
which can be derived from it should have a relativistic supersymmetry stemming from the 
covariance of the original formalism. It would be very nice if this happens because it would 
indicate that the relativistic supersymmetry is a universal phenomenon. Moreover, because 
of the geometrical meaning of the Qh mentioned before, it would also indicate that the 
supersymmetry is present because of some deep geometrical reason. This reason is because 
we basically wanted to reformulate the original field theory in a coordinate-free form, like 
the Lie derivative is. This may also explain why the partner Dirac fermions, which in the de 
Donder-Weyl formalism may turn out to be the differential forms T of the CPI for bosonic 
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field theories, are not fundamental fields but objects derived from the original fields of the 
theory. This may also explain why we do not see them in nature as independent degrees of 
freedom. The same trick would apply if we start from a fermionic theory and do its associated 
CPI. The forms would be bosonic and we would not see them because they are just differential 
forms. 

7 Conclusions 

In this paper we have extended the CPI formalism from the point particle p] to a gauge field 
theory. We have shown in details that, like in the quantum case, a gauge-fixing and a FP 
determinant are needed also in the CPI . We have also proven that this approach to the CPI 
is equivalent to a different one built out of the BFV formalism. In this first paper we have 
concentrated in building rigorously the formalism and the various universal charges that are 
present. In later works ^^'CHl we will apply the formalism in the directions outlined in the 
Introduction and in Sec. 6. We have preferred to divide the project in this way because the 
formalism of the CPI for YM is rather heavy and we could have run easily into problems if we 
had tried to do too many things without having laid the foundations properly. Last, but not 
least, the reader may wonder of what happens in the CPI to the "minimal coupling" scheme 
which was needed at the quantum level to couple particle degrees of freedom to gauge fields. 
This issue has already been thoroughly worked out in Ref. j32l to which we refer the reader 
for details. 
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Appendices 



A Construction of the canonical action of Eq. (\2A\ ) 



In this Appendix we want to derive in detail the "canonical action" (|2.4|) from the Lagrangian 
density C = —-^Fj/I^^F^jj. First of all, from C we construct the energy-momentum tensor: 

dC 

The Hamiltonian is the integral over the space variables of the component G^'^ of such tensor, 
which is given by: 

000 ^ F^,\d^Al) + \FrF'^^ 

= vr, • > - ivfa • vr" + ^Ba ■ B\ (Al) 

where we have used the definition of the conjugate momentum vr^ = —F^^ and of the "mag- 
netic field" = -e'^^Ft. So we get the following expression for the Hamiltonian: 

H, = j dxeOO = y dx {^TTa ■ > - ]fa ' + ^^a " • (A2) 

Taking into account the definition of the momentum: 

< = -FT = K + + C^-AlA'^, (A3) 
we can easily derive the equations of motion for Aa- 

L = ^a- + C^'A^Al. (A4) 
Using the previous equation into HA2|) we get: 

+2- 



The canonical Hamiltonian He is not uniquely determined because of the constraints derived 
from the gauge invariance of the theory. First of all, from (|A3|) we have that the antisymmetry 
of the tensor Fa'^ and the definition itself of conjugate momenta imply that 

7r0 = 0, a = l,...,n. (A5) 

This means that we have a number of primary constraints, equal to the number of generators 
of the Lie algebra that we are considering. It is clear that we can add to He an arbitrary 
linear combination of the constraints and consider the following function: 



H' = y dx {He + 7;?7r°) 
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as generator of the time evolution. The equation of evolution of the primary constraints 
produces the following relation: 



(A6) 



From the consistency condition, i.e., the requirement that the primary constraints are con- 
served in time, we get that the RHS of ()A6|) must be put equal to zero: 



(A7) 



The n secondary constraints (|A7|) . together with the n primary constraints ()A5|) . which we 
have obtained before, form a set of 2n first class constraints.* An arbitrary combination of 
the secondary constraints can be added to H' to get the following total Hamiltonian |,18) : 



(A8) 



= I dx Qvr, • Tf" + ■ B- + <7r0 + {A^'^ + vt){V ■ + C\,^, ■ A')^ . 

Starting from H-j. we can derive the equations of motion for the field A(x): 
' Al{x) = {Al{x),R^] = -v^,{x) 

< 

, Ab = ^b + C^l (AO" + ) 1, - V (AO + v^b) ■ 
Comparing the result with ()A4|) we can put f2a=0. Replacing this result into ()A8|) we get: 

from which we derive the canonical action of Eq. 1)2. 4|) 



dx ( ^TTa • Tf^ + \^Ba ■ B'^ - i°vr°'' + ( V • - C^VTfe • A< 



dx 



X • vf" - ivr^ • vf" - ■B'^-Al (-V • fr, + C^,^, ■ A' 



First class" basically means that the Poisson brackets of two constraints produce a linear combination 
of the constraints themselves, so this combination is zero on the constraints surface. 
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B Proof of Eq. (EUni) 



As we have seen in Sec. 3, in the case of YM theories we deal with mixed variables, so we 
have to consider the superdeterminant of the Jacobian of the transformation that brings us 
from the solutions of the equations of motion to the equations of motion themselves. How 
can we exponentiate such a superdeterminant? In the case of purely Grassmannian even 
variables one considers the Jacobian of the transformation, which can be exponentiated via 
a pair of Grassmannian odd variables: 



det" 



V-TTaV^b exp 



In the case of purely Grassmannian odd variables one considers instead the inverse of the 
Jacobian of the transformation, which can be exponentiated via a pair of Grassmannian even 
variables according to the following rule: 



det M' 



ab 



VOaVOk exp 



-OaM'^^'e^, 



It is quite natural to think about the possibility of exponentiating the superdeterminant via 
a generalized Gaussian integral. In particular, we expect to have a number of Grassmannian 
even auxiliary variables which is twice the number of the Grassmannian odd equations of 
motion and a number of Grassmannian odd auxiliary variables which is twice the number of 
the Grassmannian even equations of motion. In this Appendix we want to prove the previous 
statement, i.e: 



sdet 



A B 
C D 



d^d^dxdx exp 



A B 
C D 



(Bl) 



where 6 and 9 are Grassmannian odd variables while x and x are Grassmannian even variables. 
Let us define the following quantities: 

a = 6^ AO, b = e^Bx, c = x^C0, d = x^Dx, 

which are Grassmannian even, since the blocks B and C are Grassmannian odd, while A e 
D are even. At this point we can expand the exponential of Eq. ()B1|) as follows: 



exp [—(a + b + c + d)] 



i-y 



^ ml 



-{a + b + c + dy 



and use the "multinomial formula": 



(a + b + c + dy 



E 



ml 



(m - n)!(ii - 12^(12 - isV-hl 



Because of the presence of Grassmannian odd variables, a lot of the previous terms turn 
out to be zero. Let us indicate with n the number of variables 9 which are needed for the 
exponentiation. For n = 2 the only terms which are non zero and give a contribution to the 
integral (|BH1 are a^, abc, b'^c^, for n = 3 the only terms which survive are a^, a^bc, ab'^c^, 
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b^c^ and so on. On the other side d can have an arbitrary exponent, since it is made up by 
purely Grassmannian even terms, so we can write: 



/ = 



dOdOdxdx exp [— (a + 5 + c + 



d6'd6'dxdx<^ e""* + (-a + 6c)e" 



+ 



+ 



1 /2! 



2! V2! 
1 /3! 



1 



3! 



3! V 1!1!1! 



abc + 



1 / 4! 
4! V2!2! 



1,2 2 
C 



1 / 4! 



3! V3!" 7^4! V2!l!l!' 



1 / 5! 



5! V2!2! 



1 / 6! 



6! V3!3! 



b^c' 



When we perform the Grassamnnian integrals, the only non-zero contribution comes from 
the terms which contain n variables 6 and n variables 6. Such terms include: 



(-) 



2n-l 



(n + 1)! Wn-l)!l!l! 
(2n- 1)! 



a'^-'bc + . 



(2n- 1)! \ l!(n- l)!(n- 1)! 



(2n)! 
(2n)! I n!n! 



(-r 



n! 



n! 



+(-) 



n-l 



(n-1)! 
n! 



(ra-2)!2!2! 



+ ... 



(n-l)!(n-l)! 
Prom the definitions of a, 6 and c we get: 



!t7.! 



n!n 



d^d^dxdx- 



+ (_)"-! 



n! 



(n-l)! 



A9)''-'{e'Bxx'C9) + --- 



-j(e^A9WBxx*cer-^ + (-r^ie^Bxx^cer 

(n — l)\{n — 1)1 n!n! 



t:*Dx 



Now we have to evaluate all the integrals. First of all, let us consider: 

{9 Bxx CO) = Oi^Bi^j^Xj^XnniCm-iOi^oi ' ' ' ^i„Binj„Xj^^Xrn„CranO„Go„- 

Since B, C, 9 and 9 arc all Grassmannian odd, every block which appears in the previous 
expression is Grassmannian even. This means that we can first of all fix the indices of the 
variables 9 and then sum over all the permutations (even) of n Grassmannian even blocks: 

(^9 Bxx C9) = nWlBlj-^Xj^Xjji^CjjiiOl^Ol ' ' ' ^riBnin-^jn-^mn^mnOn^On- 

Then, if we want to put in ascending order also the indices of 9, we have to change some of 
the anticommuting parameters and give the right sign to each permutation: 



{9*Bxx'C9Y = nWi9i ■ ■ ■ 9n9n |(t|(Bxx*C)i^(i) • • • (Bxx*(7)„^(„ 
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Let us suppose we have brought D to a diagonal form (we can always do it via a unitary 
transformation), and consider the bosonic integral: 



(B2) 



Let us perform the following change of variables: 

The integration measure in HB2|) becomes: Jj2dtjdwj, while the exponential takes the form: 



exp 



. These expressions are both symmetric for the exchange t ^ v. This 
means that, to give a non zero contribution to the integral, the coefficients of the exponential 
in (|B2|) must also be symmetric for this exchange which is equivalent to the following one 
between the original variables: 



From these considerations we deduce that only the terms of the type: 



XiXi) ' exp 



can give a non zero contribution to the integral ()B2|) . The associated integrals can be solved 
in the following way: 



dxjdxj JJ(xiXi)"' exp 



Y^x'^Dkkx'' 



n 



rii 



dDiiJ \Di 

Let us notice that if D is a diagonal matrix then the following sum: 

^in A^jiCiia{i)-B2i2 A^]2Ci2«T(2) • • • det(L'"^) 
contains terms which are proportional to: 

n9idet(D-^), 

where rij is the number of eigenvalues Da of the matrix D. This is just the result of the 
integration HB3|1 . So we can write: 

J dxdxie'Bxx'COYe-''"^'' = n\det{D-^){¥BD-^CdY . 
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Similarly we have: 



Finally, collecting all these results, our integral is given by: 



I 



dede- 



(n-1)! 



{e'AeY-\e'BD-'ce) + 



— — 77 ' — 

-(e^AeWBD-^cer-^ + (-r—(9^BD-^cer 

[n — 1)! n! 
^J-{9\A - BD-^C)eTdei{D-^). 



dOdO 



det{D-^) 



Remembering the standard rule of exponentiation of the determinant we get: 

/ = det(^ - BD-^C)det{D-^). 

The last expression is just the definition of the superdeterminant. This completes our proof 
ofEq. (HnKl) . 
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C Equations of the Jacobi fields 



In this Appendix we would like to prove that Eqs. H3.2U() are identical to the equations of 
evolution of the first variations d£^^ of the fields i.e. to the equations of the Jacobi fields 
P]. For convenience, we will split in two sets and we will indicate the anticommuting 
variables with lower capital Latin indices and the commuting ones with Greek indices. The 
first variation of the action is: 



5S 



dx 



Modulo some integrations by parts, we get that the second variation of the action is given 
by: 



which gives the following equations for the Jacobi fields 5^^: 

r uj^,5C - 5ed^,daH - di'^d^d.H = 0, 
1 toab^i^ - 5edadf,H + 5i^dadbH = 0. 

Multiplying the first equation by lo^^ and the second one by u^"" we get the following equation 
for the first variations 5^: 

In the previous equations all the derivatives are right derivatives. Using the relations: 
we can rewrite the equations of motion for as: 

By comparing this equation with (|3.2U() we can conclude that also for YM theories the fields 
r evolve in time just like the Jacobi fields 6^. 



d^H = Hd^, 
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D Proof of Eq. (001) 

In this Appendix, making explicit use of the equations of motion, we want to prove in detail 
that L^^^ = L^P, Eq. jUni), i.e. that: 

-vr.V . + A^vr, + > . vf^ - ^vr, . vr'^ - . - AV, + TO, + C^Pa 

First of all, let us consider the gauge-fixing (GF) parts. From the LHS of (|D1|) we have that 
the GF part of L^^^ is given by: — tToV . A"" + A'^tTq. So, using our conventions 9'^ = (— V), 
this GF part of the BFV Lagrangian density can be rewritten as: 

- TTaV ■ A^ + X'^TTa = "VTaV . A^ + TTadtA'', 

= -TTadf^Al (D2) 

where we have used the fact that the Lagrange multipliers A" are the time components of the 
fields A"". What we have obtained as last term of (|D2|1 is just the GF term of the RHS of Eq. 
HD1|) . Now let us compare the terms of the Lagrangian densities, L^^^ and L^^ , associated 
with the bosonic fields. Here we can use the equations of motion of the fields A"". 

= ^« + VA" + Ca^X'^A" 
and the definition of the secondary constraints: 

aa = -V-TTa + C^^^vfft . A" 

to get 

> • vf„ - ivr, • vf" - ^Ba ■ B'' - AV„ = fr« . vr, + VA'' • tt, + C",,A'^1^ . vf, 

-ha ■ Tf-^ - ^Ba ■ ^" + A"V . na - X'^C^.^b ■ A^ . 

Up to surface terms, the previous identity becomes: 

A'^-^a- ■ vf" - \Ba ■ B'^ - AV„ = 1 (vf„ . vF" - . B'^) . (D3) 
From the definition of the generalized magnetic field B^ we have 

Bl = le'^^F'J =^ -\Ba . B'^ = -\f^F^. (D4) 
Similarly, from the definition of the conjugate momentum we get: 

< = -FT ■ = -\fTfSu (D5) 
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Collecting together HD4|) and (|D5|) we have: 

which, once it is replaced in (|D3|) . gives: 

> • - ■ vr" - ^Ba ■ - X^aa = -\Fr F^,. 

This proves that the terms associated with the bosonic variables on the LHS and on the RHS 
of Eq. ()D1|) are the same. 

Finally, we have to prove that the terms which refer to the Grassmannian odd variables of 
L^^^ and L^^ are the same. Here we need the equations of motion of C and the definition of 
the covariant derivative, i.e. 

Using them in the L^^ density we get for the Grassmannian odd part of L^^: 

-id^^CaD^.^C' = itaC'' - itaClJO'^Al - id^CaD^.^C'. (D6) 
Up to surface terms, the correspondent terms in the L^^^ of Eq. 1)3. 8() are given by: 

iCad^Di,,^C' - iPaP" + X'PaCl,C' + P'^Ca + C'^Pa 

= -id'^CaDkbC' - iPaP" + X^'PaCl.C - P'^t a + C'^Pa 

= iCaC" - iX^'taClfi' - id'^CaDj^.^CK (D7) 

Expressions (|D6|) and (|D7|) are equal. So, since we have proved that also the terms which 
refer to the Grassmannian odd fields are equal, we can conclude that on shell the BFV and 
the FP Lagrangian densities are equivalent: 

rBFV rFP 
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E Equivalence of the FP and the natural CPIs 



In this Appendix we want to show that it is possible to integrate away the extra fields present 
in Zq^i^^^ to get just Z^^'^^^' and this proves the equivalence of the two approaches. As 
we announced in the paper, we will divide the proof in three parts, showing separately the 
equivalence of the GF, the ghost and the bosonic parts. 



E.l The ghost part 



First of all, we want to compare the ghost parts. By making explicitly the Grassmann 
integrations over 9 and 9 in (|4.15j) we get the following dependence of on the fields of 
the theory: 



t^'' Pb 



(El) 



Going to the functional integral, we have that its ghost part is given by: 

VCVKpVCVKpVT^VT^VTpVTp exp 
For the "natural" CPI of Eq. (|2.12l) the ghost part is instead given by: 



;7YM(FP) 



I I dx£gh 



,7YM{NAT) 
^gh 



D^Dt] exp 



(E2) 



(E3) 



We already know from (|4.7j) that ip can be identified with C; if we identify Ap with rj, we 
get: 



J Vi/jVt] exp j dxd^i)aD^i^i]h 



VCVA-pexp 



i I dxid'^CaDll'Ap^ 



(E4) 



It is clear that in the path integral Z^'^^' of HE2|) there are six other functional integrals, 
in the variables C, Ap, F'^, F'-^, Fp, Fp, which are not present in the associated "natural" 
path integral of Eq. (|E3|) . In what follows we want to show that, performing explicitly these 
integrations, we get a quantity which is independent of the fields of the theory. Let us start 
by evaluating the integral (|E4() over Ap: 

T = j VAp exp jy da; {-d^" Apad^.C' + ApaC^^^^A^C^)^ . 

After an integration by parts we get: 

T = I VAp exp I [ dx \Ap^ (d^'d^C' - C^d^" {A^^C^) 



~ 6 {d^d^St - C^d^A^ - C^^^Ald^C' . 
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Using the gauge-fixing condition d^A^ = and the properties of the Dirac deltas, we can 
write 

T ~ 6{C'') detid^'df.S^ - C%A''^^^'). (E5) 

The Dirac delta 5{C^) allows us to perform immediately the integral over C, putting C = 
in all the terms appearing in (jEll) . Next we have only four other functional integrations to 
perform: 

detidf'd^d^ - C^AI^df") j VT^VT^VTpVTp exp j dx (id^'C aC%T^'^^Tp^ 



Using the definition of covariant derivative and integrating by parts, we can rewrite one of 
the terms appearing in the weight of HE6|) as follows: 



Note that also the integral over r*^° produces a Dirac delta: 

6 [{d'^d,6i - c%A<;^dnrp^] = det-i {d'^d.dt - c%A';^d^^) s{Tpj. (et) 

The presence of the inverse of a determinant is due to the fact that the CPI variables F 
associated with the BFV ghosts are bosonic. The Dirac delta S{T-pJ present in HE7|) allows 
us to perform the integral over Fp^. The result is that we can put everywhere: F73- = in 



HE6|) . The two determinants of Eqs. HE5|) - HE7|) cancel each other and the result is: 

I VT^VTp exp |i I dx [d'^CaC^T^.T^' - d'^Tpad.T^" + S'^Fp^C^^.A^F^") | 
= det-\d^d'^6t-C%A';^dn. (E8) 

So the only difference between Z^'^^' and ^J^t^*"^' ig given by the determinant in (|E8|) . 
We want to prove that such determinant does not depend on the fields A'^. A functional 
determinant must be considered as derived from a matrix which depends on the space-time 
variables x: 

D = det-^ {d^d^'S^Six - x') - C%6{x - x')A''^d^) . 
Using the relation 5{x — x') = d^9{x — x'), we get: 

D = det-H5^)det-^ [5^5{x - x') - C^bOix - x')A''^] det-^(a^). 

The terms det~^(9^) and det~^{d^) are independent of the fields A. So it remains only to 
prove the independence of the term: 

det-i [6^6{x - x') - C%9{x - x')A1^] = exp-^ Tr ln(l - A) 

= exp-^ Tr [-A - A'^/2 ■■■], (E9) 
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where A = C°^^6{x — x')A'^. On the RHS of HE9|) we have only performed standard manipula- 
tions. The trace in ()E9|1 is zero, thanks to the antisymmetry of the structure constants C^^. 
Consequently, the determinant appearing in ()E8|) is independent of the fields of the theory. 
This completes the proof that the ghost terms of the two functionals ^ym(nat) ^'ymcfp) 
are equivalent. 

E.2 The bosonic parts 

The equivalence of the bosonic parts is much easier to prove. In fact from Eq. (|2.14|) we have 

The previous Lagrangian was obtained by exponentiating the equations of motion derived 
from the canonical Lagrangian (|2.4)) : 

L = 1, • 7f" - ^TTa • ^» - \fI^F^^ - X'^i-V ■ na + C^.^b ■ A'^). (ElO) 

Using in HE10|) the equations of motion: 

A" = 71" + VA" + C^A'^l'^ 

we get 

i.e., just the original Lagrangian for YM theories. Using this result and the usual connection 
with the superfields we can write, up to surface terms 

Note that the RHS of this equation is also the bosonic part we would get in Zq^i^^\ So this 
proves that the bosonic parts of Z^^i^^^ and Z^^^^'^'' are the same. 

E.3 The gauge-fixing parts 

Finally, we have to prove that the two gauge- fixing parts are equivalent. In C^^'^ we have 

only one term: —iiad^A'^, while, using the standard superfields, the gauge-fixing part of C^^ 
is given by: 

= i j dede[-^ad^Ai] 

= ^ad^Ki + ^r^-a'^r^M - iTxa^^'T^t^ + kxc^^^'Al. (eii) 

We can now formally identify the fields tt^, which appear in >Cgf*^, with the fields A^o, which 
appear in HE11|) . Then, integrating Z™'^^' over Tr^, T'^" and T\a^ we get a term proportional 
to: 

6{d^K.)6{d^^T^.)5{d^T^^ = detid^6t)S{A^.)6{T^.)5{r^h (E12) 
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The reason why we obtain a determinant and not an inverse of a determinant in HE12|) is 
because there are two anticommuting fields F^m and F^m for each bosonic field A^m. The 



deltas which appear in (|E12|) are not totally painless to work out since they contain fields 
which enter the expression of /^bos- particular, these Dirac deltas make C^^^ = 0. This 
could have been foreseen: in fact, our approach is on-shell and, using the equations of motion, 
we have: 



This tells us that the Dirac deltas of Eq. ()E12|) are somehow superfluous. For example, let 
us consider the Dirac delta over F'^''. In the bosonic part of the natural Yang-Mills CPI 

•YM( 

bos 



^YM(NAT) ^j^^ terms depending on F'^m are the following ones: 



/ 



ppA^Pp^ exp 



,_d d_ 



H F^- 



(E13) 



If we call F'^'' = y then the integral in HE13|) can be written in a compact form as J dy 6[f{y)], 

with f{y) = at y = 0. As we have seen in HE12|) . the GF part produces another Dirac delta 
S{y). This Dirac delta is superfluous if the following equation holds: 



dy5[fiyMy)r^ Uy5[f{y)] 



(E14) 



The "~" sign means that the only difference between the LHS and the RHS of ()E14|) is given 
by terms independent of the fields. 

Using the following properties of the Dirac deltas: 



5{y - a)5{y - b) = 6{y - a)6{a - b), 



dy5[fiy)] 



dy 







dy 


y=0 



we can write: 



I dy6[f{y)]6{y) = J dy 



Siy)d{y) 



df{y) 




dy 


y=0 



m / dy 



df{y) 




dy 


y=0 



5(0) / dy6[fiy)]. 



So 5[f{y)]6{y) and S[f{y)] differ by a 6{0), which is independent of the gauge fields. This 
concludes our proof. 
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F Proof of dSZHD 

In this Appendix we want to prove that if two functions Oi and O2 have zero Poisson brackets 
in the standard phase space then the associated 



(Fl) 



commute in the extended Hilbert space underlying the CPI. More generally, we will prove 
that if {Ol, 02}pb = O3, then, using the commutators of the CPI and the definitions HF1|). we 
have [01,02] = iOg. 

Since each Oi is the sum of two different terms, the commutator between Oi and O2 will 
be given by the sum of four different commutators. Let us analyze them one by one. The 
one between the first terms is: 



-iAoUj''^ dj^d^O^d ^0^^Afl(_)SO2+OiO2+Ai3+A^ 



(F2) 



where each letter in the exponent of the minus sign indicates the Grassmannian parity that we 
have to introduce to take into account the fact that the CPI fields have mixed Grassmannian 
parity. The term on the RHS of ()F2|) must be compared with the Liouvillian associated with 
the standard Poisson brackets between Oi and O2: 



A^uj^^d, 



Ol d c^^° d 0O2 



A^cj^^ d sOi d ouj''^ d BO2 + A^w-^^Oi d cuj^° d b d oO 



\{Oi+C)B 



iA^Lo^^ d sOi d o^^^ d bO, - iA^u""^ d^d^O.d ^OiW-^^(-)^°2+°i''2+^^+^ ) , 



which is just —i times the RHS of ()F2|) . Similar relations hold also for the terms coming from 
the other commutators, i.e.: 



Aauj^^ d flOi, iFc 0;°^ d BO2 d 



\02 



= r^w^^ dj^dpO^d ^r^w^^ d ^o^(_)^o+02+02fl+020i . 

iF^cj^^^sOi a"cr^, Ao iv'^^'dEO^ 

= -r^a;^^?^Oiac9^^^^5^02r^(-)°^<^+^'+°^ 

+r^u;°^6's02 5"^u;-"^6'sOi9"cr^(-)°20i+Oi+02. 



Let us consider the second term coming from the definition of O3, i.e. 



Ol 9 c uj^"" d 



\0i+02 



(F3) 
(F4) 
(F5) 
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From this equation we get four different terms because the derivatives w.r.t. and can 
be apphed on both Oi and O2. Let us write exphcitly these four terms with the correct 
grading factors. By applying d b over Oi and d e over O2 we get: 



01+02 



which, modulo a factor i, is just the first term on the RHS of HF5|) . By applying the derivatives 
w.r.t. and over O2 we get 



which is just the RHS of Eq. ()F3|) multiplied by —i. If we apply both the derivatives over 
Oi we get: 

which is —i times the RHS of Eq. (|F4|) . Finally, if we apply the derivative w.r.t. over Oi 
and the derivative w.r.t. over O2 we get: 



which is —i times the second term on the RHS of ()F5|) . In this way we have proved that if 



{Oi,02}pb = O3 then the associated Lie derivatives Oi and O2 satisfy 



Oi,02 
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G Proof that 



= 



In this section we want to prove that the two BRS charges, the universal Q of the CPI and 
the one Vt, due to the gauge invariance of the theory, commute. First of all we need the 
explicit expression of the BRS-BFV charge, which from Eq. ()5.8|) . is: 







1 ._p 

p, + 2- 
1- 



The expression of the CPI BRS charge is much more simpler (|5.ip : 



Q 



ir^A^ = ir^fcA^" + «r^° A^fc + ir^ Aa. + ^r^^A^ 
+zr'^"Ac<- + ir^^Apa + ir^"A7=. + ir-^^Ap . 



(Gl) 



The non zero commutators of the theory are: 

r,A^]=.<5f, [^^^^] =5f. 

Let us start by rewriting the non zero commutators of the first four terms of with Q: 



d,T^^ A 



(Ofc^^.)^^^^^^?A 



(9feAA^)C^^^^ Act 



{duAAi)T<'\ 



The first and the second commutator are opposite, as well as the third and the fourth. 
Something similar happens for all the other commutators. A first set is the following: 



-C\AAlAp^,iT<A^,^+iT'''Aj,, 
t&^A^^lTp^^iVP^Ap^+iT^'^A^. 



fib _ipAjA_ 



~tb k 
' ac^b 



-c"LX'A,jC^ir^ Ac. + ir-dA^, 
iCtr-^'^^,rp^,ir^'*Ap^ + zr^.'A^. 



ac 



iC\,Y<Y^kC\iT<A^,^iV^'Ac. 
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In the previous formula each commutator appears together with its opposite, so they do not 
give any contribution to the commutator Let us go on with other commutators: 



-ctr^^^c'^A^. -ictr^.r^^r^", (G2) 

k k 

-Ci,A^.AlT^"' + Ci,A^^T^lC^. 



Also this block of three commutators gives a zero contribution to the commutator [Tl, Q] since 
aU the terms in the RHSs of Eq. are one the opposite of the other. Let us consider the 

next block of four commutators: 



-iP"AAa,ir^ Apb 
T-^J''\iT^''A-^^ 



ir^" Aa. , [-l^Ta. , ir^'AJ = -ir^" Aa. 

-zr-A^^, [zA^^vr„,ir-^A.J =zr-A^^. 
Finally, we have the last block of three blocks of three commutators: 



1.^ 
2' 



Pc 



-iCl,rP^C'Tp^,iTP^A-p^ + zr^-^A^, 



1^ 
'2 



=^PaCl,C'Ap^,ir'''^Ap^ + iV'^'Ac. 



1 



=-C'lTP-C''Ap + ^Cl^PaT'^"Ap 



1 ^ 



a -3 T-C" 



1 — - 



--iPaClj:pT^\iT-^Ap^ 



-irP^C%j:pC^,iTP^Ap^+iT^' Ac. 



1^ 



;-PaC%,ApC', iT^'^Ap^ + iV^'^Ac. 



be'- 



1 ^ 



-c",,r^«Ap^c'= + -c-,,p,Ap^r 



^iFc'^Cl.C'T^^iT^'Acd 
-UTc^C'l,T^''C'^,iT^'Acd 



--Ac^Cl,C''C'^,iT^''Acd 



-C 1 Ac« - -iC f,^l c^L i 

-Iac-^cit^'c^ + Iac^cit^^cK 



It is easy to realize that the RHSs of each of the three blocks of commutators contain terms 
which are one the opposite of the other. This means that the contributions to the commutator 



[fi, Q] of each of the three blocks above are zero. This concludes our proof that Q, 



0. 
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